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Abstract 

The billiard problem of statistical physics is considered in a new geometric approach with a 
symmetric phase space. The structure and topological features of typical billiard phase portrait are 
defined. The connection between geometric, dynamic and statistic properties of smooth billiard is 
established. Other directions of the theory on development are pointed out. 

1. Introduction 

Billiard is one of the most important models of statistical physics and chaotic dynamics. G.D. 
Birkhov suggested regarding billiard as a typical conservative system [1]. A.N. Krylov based his 
explanation of solid spheres gas statistic properties on exponential divergence of its "billiard" 
trajectories [2]. In the works by Ya.G. Sinai and L.A. Bunimovich on phase trajectories mixing 
in scattering and defocusing billiards Boltzmann's hypothesis of molecular chaos found its further 
grounding [3] . Now billiard became a paradigm of deterministic chaos [4] of classical systems and 
is often applied [5] for the research of their quantum "twins" . A lot of applied physics problems 
can be reduced to a billiard problem [6]. 

A classical billiard problem is in studying its character and distribution of its trajectories. 
Among the typical billiard motions one can point out the following: periodical, quasiperiodical 
(integrable) and irregular (chaotic) motions. Compound billiard dynamics appears in the phase 
portrait structure of the corresponding mapping. The latter is plotted using different geometric 
methods or Poincare sections. For the specification of a billiard ray it's usual to choose local 
Birkhov coordinates: natural parameter I in the reflection point on the border of the billiard d fl 
and the incidence angle 9 in the same point. They stand for canonical variables - the coordinate and 
the moment for Hamiltonian description of the system. Many important properties at this choice 
of phase space coordinates stay unnoticed. Let us choose another unifying approach. It identifies 
billiards with reversible mapping (with projective involution) in a symmetrical phase space. In 
its framework one can join together geometric, dynamic and statistic properties of billiards. Such 
fundamental mechanisms are analyzed in this paper. 

2. Symmetric Coordinates 

Let us describe geometric propagation of the rays of billiard (Fig.) as a reversible mapping 
B of the phase space Zwith symmetric coordinates (zi,z^). The pair of these coordinates defines 
two successive reflections of a billiard ray from d f2. At the same time, each of the coordinates 
corresponds to some parameterization of the billiard border, r\ dn = r(z) — (x(z),y(z)). The 
following topological construction appears: Z oc d O X d ft. For a closed planar billiard one can 
accept z G S 1 (circle) or z € I = [0, 1] the periodicity being r(z) = r(z + 1). So we'll have a phase 
space as a torus Z = T 2 = S 1 x S 1 or its unfolding II = I x I on the plane. After each reflection 
of an arbitrary (incoming) billiard ray with the coordinates (zi, Z2), we have a (reflected) ray with 
new coordinates (z[, z' 2 ). As a result, the evolution of these successive reflections is described with 
a billiard cascade {z\, Z2) —> (z[, z' 2 ) of the form [7] 
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z[ = z 2 

4 = f{zl,Z2) 



; f(f(zi,z 2 ),z 2 ) = zi;g(f,z 2 ) = g(R(zi, z 2 ), z 2 )); R(z,z') 



a(z r )z + b{z') 
b(z')z - a(z') 
(1) 



with the involution / = f(zi,z 2 ) (on the first argument^), which is defined by implicit depen- 
dence on the corresponding fractional rational involution R, R(R(z, z'), z') = z. The coefficients 
a(z) = n^(z) — Tiy(z); b(z) — 2r? x {z)r? y {z) are expressed with (Cartesian) components of the 
exterior normal field n(z) = n ext \dn on the border d fi, ft = (n x ;n y ) = (y'(z);—x'(z)) (the 
stroke marks differentiation). Function g depends on the form of d £1 and g{z\, z 2 ) — g(z 2 , Zx) = 
— [x(zi) — x{z 2 )]/[y(zi) — y{z 2 )\. The choice of Cartesian coordinates is rather relative, because 
every covariant substitution Z\ — > ^B(zi); z 2 — > £B(z 2 ) preserves the form of the mapping (|l|). The 
mapping built can be used for the description of billiard with a border of the most general type. 
During its derivation only the condition of elastic reflection and no restrictions of smoothness, 
curvature, convexity, simple connectivity of the border and so on were introduced. 

The mapping (|l|) is invariant to the substitution, S := z\ — ► z 2 \ z 2 — ► Zi, of the incoming 
ray to the reflected one, i.e. B o S = S o B for the composition of transformations. This means 
reversibility of the constructed mappings. The physical reason of this is reversibility of the system 
to the changes of the time sign (the direction of the motion). This is a global property. In 
the billiard cascade, phase trajectories with opposite directions of the motion or with opposite- 
directional initial rays (210,220) e (z 2 q, zio)are present simultaneously. This requirement of local 
reversibility is stronger. The inverse of the ray reflected with its successive reflection makes the 
initial incoming ray. Mathematically it leads to the appearance of a involution / in the mapping 
(0) . The symmetry (reversibility) leads to the symmetry of the phase space and the phase portrait 
of the mappings (|l|) . For every element Z there is one symmetrical to it relative to the diagonal 
A = {(zi,z 2 ) £ Z | z\ — z 2 }. For every function ^onZ 



That's why it's natural to regard the coordinates of Z as symmetrical. This symmetricity essen- 
tially sets them apart from the variables of Hamiltionian description of the billiard. In particular, 
the locality of Birkhov coordinates causes explicit mappings to be obtained only for the simplest 
geometry d £1. In non-local coordinates, (zi, z 2 ) this problem can have a general solution. 

In the research of periodical trajectories of the billiard the powers of billiard mapping B k are 
also used 

B k := {z[ = fk-i(z,z);z' 2 = f k (zi,z 2 )}; fk(zi, z 2 ) = f(fk- 2 (zi, z 2 ), fk-i(zi, z 2 )) . (3) 

They include billiard "compositions" /j., where k = 0, 1, 2, . . . ; fi = /; /q = z 2 ; f_x = Z\. They 
lose the property of involution, but preserve reductibility to fractional rational transformations. 
The mappings (||) describe "pruned" billiard trajectories with the omission of a set of (k — 1) links 
(successive reflections) . 



All the geometric properties of a billiard are established in the specialization of mappings (Q) . 
They are concretized in the features of the involution/. In the appropriate (local) coordinates it 
can be reduced to fractional rational involution R. Projective transformations are described with 
fractional rational functions. Billiard is one of those transformations. In every reflection point 
incoming and reflected rays are joined together by a harmonic transformation G. For G projective 
invariant (a complex relation of four rays , incoming i, reflected r, normal n and tangent t) is equal 
to (i,r,n,t) = — 1. In geometric terms involution looks the simplest 



X{zi,z 2 ) = x{z 2l z x ) 



(2) 



3. Billiard Geometry: Involution Properties 



G(i;n,t); GoG = Id 



(4) 
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Id is an identical transformation. Let us emphasize the locality of the projective property 
of the billiard. The concrete form of G depends on the guiding-lines of the normal in the point 
of reflection, that is, on the form d f2. Harmonic mapping is an involution and changes the 
sequence order of the ordered projective elements to the opposite. The monotony of / on the first 
argument is the consequence of this. This monotony of piece- wise continuous / (involution can 
stand discontinuity) is true for every billiard. Using the correlations (|l|), the involution can be 
laced of local branches of the form / = g^ 1 o Ro g (the choice of <? _1 branch is dictated by the 
variational principle, the minimality of distance between the points of reflection). From that the 
property of monotony immediately follows 



df(z 1 ,z 2 )/dz 1 <0 . (5) 

Fractional rational functions are dense everywhere in the space of continuous functions. In 
fact, this means the possibility of arbitrary precise approximation of different physical systems 
with their billiard models. This fact is used, for instance, in the analysis of energetic spectra of 
multi-particle systems, the description of kinetic properties of continuum (Lorenz gas model) etc. 
If sinus and cosine have physically appeared from the problem of oscillator, then fractional rational 
functions can be generated by billiard. 

The reflection of ray beams from the border of the billiard can be of diffractive, focusing and 
neutral character. This depends on the curvature of d fL The representation (|I|) gives the following 
property 



^§^}=sign{k(z 2 )} , (6) 

where K is oriented curvature in the point of reflection. For the convex border K > involution 
appears to be a monotonous function on both arguments. (For instance, for a circle, f(zi, z 2 ) = 
2z2 — z\ (mod 1).) On a torus, Z = T 2 , such involution has no breaks (Unlike everywhere 
dispersive billiard, K < 0, with lacunas in phase space.) 

Involutivity and projectivity are the main geometric properties of a billiard. The geometry 
(form) of its border defines the explicit form of involution. At the same time, it also defines the 
dynamics of the billiard. 

4. Billiard Dynamics: the Structure of Symmetric Phase Space 

Let us analyze the structure of symmetric phase space of a typical billiard (see Fig.). This 
principally solves the question of the types of dynamics and stability. For more simplicity let us 
regard the border d fl as a differentiable as many times as needed. The research of billiards in 
polygons (without curvature) needs some peculiarities. 

For high-quality research of phase portrait of the mappings and its local bifurcations normal 
Poincarc forms are especially useful [8] . In the symmetric approach the theory of normal billiard 
forms appears to be the most advanced. This is connected with the flexibility (a wider class of 
allowable variables) of reversible systems. Any changes of variables in Hamiltonian approach are to 
preserve the conservation character of the mapping with the Jacobian J = 1 (canonical changes). 
Whereas the mapping ([!]) doesn't demand it. It Jacobian J = — df{z\, Z2)/dz\ > can take 
arbitrary values, < J < 1; J > 1. As a weak limitation, the demand for the mapping (^) to 
preserve measure remains. This means that J — J(z) — p(z) / p(B(z)), where z = {z\,Z2)\ Bz — 
(z2, f(zi, Z2)) should be true. The proof uses the equation of Frobenius-Perron for the density p of 
invariant measure (see further) and the symmetry (^) for it, p{zi,Z2) — p(z2,zi). This limitation 
can always be met preserving the main property of involution / o / = id in new coordinates. 

Omitting the details, let us present the expression for universal normal billiard form in symmet- 
ric coordinates. It is true in the neighbourhood of an arbitrary cycle of porder (periodic trajectory 
of ^period) 

NB P = { Z ' 1 = ~^p~ iZi + ^p-i^ + z 1 P(z 1 z 2 ) 

(4 = -Mp^l + V P Z 2 + Z 2 Q{ziZ2) 



sign 
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where the coefficients of the linear part are defined by the expansion of "compositions" and 
fp (see formula (0)) in the initial point neighbourhood of the cycle under consideration. They 
constitute the matrix of L linear part. Its determinant is equal to one, det L = 1, that agrees 
with conservation of measure and the condition J(C) = 1 for every cycle C of mapping (^). 
Homogeneous polynomials P, Q (without absolute terms) define nonlinear additives. Their explicit 
form depends on the involution of billiard /, that is, on the form of d fl. 



The character of the cycle depends on the size of trace trL. For an elliptic cycle 



tr L 



< 2, 



for a hyperbolic one tr L > 2. In the neutral case, for instance, for a billiard in a circle, trL = 2. 
It can be shown that for any cycle, corresponding to a periodic trajectory, passing through a 
concave section with concavity Kiz^) < 0, trL < —2 will be true. That's why the trajectories near 
such cycles always are unstable and exponentially diverge from one another. Near elliptic cycles, 
including 2-cycles, regions of regular motion form. With the loss of ellipticity they are ruined, 
first forming stochastic layers and then, when the latter are covered, a chaotic sea. In reality 
this mechanism looks much more complicated, for example, with an intermediate pass through 
Cantor-tori etc. Normal forms (0) let us trace typical properties of such bifurcations, taking place 
when the billiard border is deformed. 

The diagonal A(zi = Z2) contains all fixed points of the billiard, BA = A. This follows from 
the diagonal property f(z, z) = z of billiard involution, resulting from its coordinate expression (Q). 
For a convex billiard in the neighbourhood of the phase space diagonal, normal form (0) can be 
reduced to the mapping of a turn, i.e. a particular case of a billiard in a circle. Here the structure 
of elliptic and hyperbolic cycles of arbitrary high order is shown. The motion stays regular. The 
appearing of negative curvature ruins this situation. There is no unified transformation (or the 
integral of motion) near the diagonal because of appearing breaks of billiard involution. 

Analytic research of the symmetric phase space structure can be continued using geometric 
methods. Here the important advantages of the new approach are seen. In addition to regular 
and chaotic components of motion the phase portrait can contain regions of forbidden motion - 
"lacunas" Land regions of degenerated motion - "discriminants" D. 

Lacunas (Fig.) appear in the billiards with regions of negative curvature. They occupy the 
phase space part, the points of which correspond to the rays lying outside of the billiard region 
fi. The coordinates of these rays meet the condition r(zi) — r(z2) ^ fi. This condition defines the 
inner region of lacuna L in Z. The form of the lacuna is defined by its border d L :— {(zi, z^) G 
Z\z\ = \{z2)}- (Another parameterization is possible Z2 = A (zj,). In this case functions \(z) 
and X~ 1 (z) specify the same simple closed curve, but passable in different directions. When one 
its branch is higher than Aand the other is lower, and vice versa.) The border dL comes to the 
diagonal A transversally and crosses it twice in the points with coordinates (zq, zq), corresponding 
to the points of inflexion, K(zo) = 0. 

The forbidden billiard rays (points of lacuna) lie in classically inaccessible region - geometric 
shade, generated by the regions K < 0. The number of lacunas (on a torus Z = T 2 ) is equal 
to the number of negative curvature components d f2. Every lacuna is a simply connected set. 
The contrary would mean non-closed character of d Vi. With the appearance of lacunas a part of 
diagonal A is cut out. The corresponding fixed points disappear. For everywhere dispersive Sinai 
billiard, lacuna absorbs all the diagonal and the mapping (|l|) will lack all fixed points. At a special 
configuration of such a border d f2 one can cut out cycles of higher order, p > 2. 

In a topological way one can glue up the lacuna on the torus with a two-dimensional manifold. 
According to the rule of d L bypass, it can be only a piece of a projective plane. This is directly 
connected with the projectivity of the billiard. On a projective plane, metric conceptions "inside" 
and "outside" of a closed region lose their sense. (For example, a closed curve and a right line 
that doesn't cross it on a plane may have common points after central projection onto the other 
plane.) That's why "forbidden" rays turn out to be involved into the general billiard flow. Such 
global motion takes place on non-oriented manifold. 

On the projective plane the initial involution / also rules the motion of the rays. Almost every 
such ray (exceptional cases are of measure null) is continued to an ordinary billiard ray, further 
dynamics of which is known. As a result of further evolution, this ray after some time will return 
to the section of negative curvature d Vt, corresponding to the lacuna under consideration. This 
is specified by the mentioned hyperbolicity of cycles that contain points on the concave border. 
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Being continued then to a classically inaccessible region (preserving the direction of motion), it 
would give a new position of the initial ray (phase point in the lacuna). A recurrent mapping 
appears. It is defines by one of the "compositions" fk, included in the equation (||), the order k 
always depending on the coordinates of the initial ray (the initial point of the lacuna). The lacuna 
plays the role of a secant for the Poincare section of the billiard flow. Similar evolution also takes 
place with other points of all lacunas. Only phase trajectories of ordinary billiard rays remain in 
this case "visible". 

The condition of "connecting" for the inner rays, that are tangent to the concave region in the 
point r(z%) and that cross d SI in the points r(z\) , r[z 2 ) outside of it, defines the border of the 
corresponding lacuna 

{f{z x )-f{z 2 ), n(z 3 ))=0; (z u z 2 )eZ\L; { Zl , z 3 ) £ L ; (z 3 ,z 2 ) edL ; K{z 3 ) < , (8) 

where (., .) is the scalar product of the vectors. Solving the equation (||) according to the theorem 
about an implicit function, we have z\ = Ai(z3); z 2 = X 2 (z 3 ). Excluding z 3 , we come to the desired 
equation z\ = \(z 2 ); A = Ai o X 2 . 

The discriminants D correspond to the zone of "stuck-together" trajectories or "non-continuable" 
trajectories that cross special (corner) points of d SI. That's why they appear in the billiards with 
straight regions of borders, K = 0. Topologically the discriminants have a lacuna-like structure. 
Their border d D is defined by 



(r( Zl ) - r{z 2 ), n{z 2 )) = ; (a*, z 2 ) E Z\L ; K{z 2 ) = . (9) 

It also can have explicit form z\ — p(z 2 ). The discriminants have regular shape - squares 
(Fig.) with a diagonal, which coincide with a part of A(zi = z 2 ) in the region corresponding to 
the straight-line component d SI. 

Lacunas and discriminants make a principal property of a symmetric phase space. In fact, they 
are filled with the rays of the billiard that fell out of its ordinary dynamics. (At their passing it's 
easy to show that the billiard involution / breaks (on the first and the second arguments), the 
breaks being different from the factor of periodicity (mod 1) and are not removed when passing 
to a torus, Z = I 2 — > Z = T 2 .) There are no such non-local elements in the phase space of 
Hamiltonian approach. At the same time, these formally hidden "topological" obstacles (Fig.) for 
the billiard flow to flow around, and the diagonal A, on which they arise, play an important role 
in the chaotic dynamics and must be included into the full description. 

5. Billiard Kinetics: Invariant Distributions 

The geometry of phase space structural elements depends on the form of the border d SI and 
(or) involution /. Let us show that in the symmetrical approach not only dynamics but also 
kinetics of the billiard is connected with these characteristics. The kinetics becomes apparent in 
the case of chaotic billiard, whose deterministic trajectories (unique groups of successive points of 
reflection) have all the properties of random sequences in the asymptotic limit of infinitely large 
number of reflections. That requires statistic description of (two-dimensional) dynamic system in 
the manner of deterministic chaos conception [9] . 

In the mixed dynamics of a typical billiard both integrable and ergodic (as a rule, with in- 
termixing) types of motion are present. In this case, the general characteristic of the trajectories 
distribution in the billiard is the invariant measure of the dynamic system. Absolutely continuous 
distributions are of the greatest physical interest. It's important for them to be individualized, i.e. 
depending on the geometry of a concrete billiard. That's why universal measures like conserved 
Liouville volume (Birkhov measure [1] for all billiards) are of little use here. A symmetric measure 
with density p (zi, z 2 ) = p (z 2 , z±) already possesses this individuality. 

From the operator equation Bp = p for an invariant measure after transformations using piece- 
wise monotony (||) we have 

p(zi,z 2 ) = p{z 2 ,zi) = p(z 2 J(zi,z 2 )) ( _ d/(zi,Z2M = p ( Z2i f(zi,z 2 ))J(zi,z 2 ) . (10) 
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In some cases (integrable billiards) one can find its explicit solution. Geometrically, p is a 
two-point density; it depends on the coordinates of two points on the border d £1. The topology 
of the direct product Z oc dfl x d£l causes one to choose a special factorized solution, p(zi, z 2 ) = 
L)(zi) (j(z2). Instead of the expression ( |To| ) we get a functional equation for one-point plane lu(z): 




to (z)=cu(f)J(z,z'): f = f(z,z') uj{z)dz = -lu (f) df , (11) 

written in total differentials. The factorization is coordinated with the symmetry of p and preserves 
its normalization ||o;|| = J Q uo(z)dz = 1. 

The physical sense of lo(z) is an asymptotic plane of billiard flow reflection points (with coor- 
dinates r(z) G dCl, z G 7). This is a truncate distribution in the sense that the dimension falls 
twice. It will be very useful in the description of physical characteristics in different billiard prob- 
lems, for instance, the "probability" of ray escaping from a fixed place of resonator, wave-guide or 
detector. Besides, it is directly connected with the involution and geometry of the billiard. After 
integrating the differential relation (|ll|) for lu we have 



Lo(z)dz = C(z 2 ); f = f( Zl ,z 2 ) , (12) 



where zq is an arbitrary initial point on d il; C(z) is the function to define. With different 
character of border d OC(z) has different forms. For everywhere convex billiard it's one can just 
use the diagonal condition f(z,z) — z, so C(z) = 2j za Lu(z')dz'. In the general case the border 
9 !] = U d£l- U dflo contains regions of positive, d 0+, negative, d and zero, d f^o, 

curvature. During the defining of C{z) the solutions in symmetric "halves" of phase space over 
and under the diagonal A, that is, in the involutionally connected regions with coordinates [z\, z 2 ) 
and (f(zi, z 2 ) 1 z 2 ) are laced. In the presence of d f2_ and d Qq components connecting takes place 
on the borders of corresponding lacunas and discriminants. Summing it up, let us set the border 
d E, that divides different symmetric components of ordinary cascade E (outside special zones) 

!Z2 , (zi,z 2 ) G A 
A(z 2 ), ( Zl ,z 2 )edL (13) 
p{z 2 ), (zi,z 2 ) G 3D 

with known dependencies in the cases of lacunas and discriminants (see above). Let us note that in 
each half of the phase space A(z) is a multi- valued function (the number of branches doesn't exceed 
the doubled number of d Q_ and d tto components, but self-intersections and multiple connection 
d E are forbidden by the uniqueness of the flow. A(z) is the functional of d f2form. Connecting on 
the border d E gives us 

/ /(A(z),*) zo\ f f(zi,z 2 ) /(A(z 2 ) : z 2 )\ 

w(z')dz' = C(z) => If - [ \ia{z)dz = Q . (14) 



\ zo AO) / \ A(z 2 ) zi 

The dependence of the initial point zq, as would be expected, falls out. The equation obtained 
lets one to restore billiard involution / on the one-point billiard distribution function wand vice 
versa. At the same time both functions are connected with the equation of border d £1 by the 
expressions (^3|) and ([[[). The billiard problem takes on a single meaning from dynamic, statistic 
and geometric points of view. 

Direct solution for won / can be obtained by differentiation of Eq. (0) 



dlnco (/) _ d2 ^ Z ^/ dzidz2 
df df{z 1: z 2 ) /d ^ df{ Zl , Z2 ) /g ^ 



(15) 
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Having performed the change Z\ — > f(z\, z 2 ) and using the relation / o / = id for calculating 
the derivatives, this equation for uj can be written as 



[in^(z 1 )]: i = --^ r + 77 i f; f = f(z 1 ,z 2 ) , (16) 

where strokes mark partial derivatives. The right part of the equation ( |l6| ) shouldn't depend on 
z<x- This limits, on the one hand, the type of involution (not every involution can be billiard, that 
corresponds to the stressed role of projective transformations), and, on the other hand, the choice 
of phase space variables, allowing factorization of two-point measure to one-point ones. Covariant 
changes, leading to "factorizing" coordinates, correspond to the choice of a certain frame and way 
of border d f2 parameterization. 

In the equations ( |l(i| ) and ( |l6| ) the densities p and ui are uniquely defined by the involution of 
billiard /. The latter is uniquely defined by the border d ^equation, according to the representation 
. The invariant measures of the billiard become its individual characteristics. In a chaotic billiard 
they acquire the character of equilibrium statistic distributions. 

So, on the whole billiard analysis in symmetrical coordinates shows that its main characteristics 
are uniquely connected with one another 



dQ, -> f(zi, z 2 ) -> p(zi, z 2 ) <-> uj (z) -» f{zi,z 2 ) — > dQ, (17) 

where arrows indicate passing from one object to another. 

One of the most designing and old problems of statistical physics is finding out the transition 
from the reversibility of deterministic motion equations to irreversibility of statistic ones, see for 
ex. [10]. Generally accepted point of view is that irreversibility appears at roughening in the 
macroscopic description of the system on the kinetic stage of evolution and is connected with the 
fundamental principle of correlations unlinking. Here usually the problem of distribution functions 
calculation with given Hamiltonian (the equations of motion) is posed. In physical applications the 
inverse problem may also appear: to restore the dynamic law for a chaotic system (not necessarily 
of mechanic origin) with known statistic characteristics. It can be of special actuality for the 
system with a small number of freedom degrees. 

Statistic irreversibility prevents the reverse of the "time arrow" , but doesn't necessarily break 
the feedback of kinetic and dynamic. A remarkable peculiarity of the billiard is the possibility to 
solve direct as well as indirect problems. The form of the border d Q defines involution, on which 
the invariant measure is calculated. And vice versa: the involution (that is, the dynamic of the 
billiard) is restored from the one-point distribution of reflections on the border. The border of 
the billiard can be restored by its involution [7]. By the way such closure is the consequence of 
geometric (projective) nature of the billiard. 

6. Summary 

In the conclusion let us note the characteristic properties of the symmetrical approach. 

1. The symmetricity of the phase space. Equality of phase coordinates rights. Non-local 
character of the geometric elements involved into the dynamics. 

2. Unification of billiards as reversible dynamic systems (mappings) with projective involution. 
Reversibility and projectivity of the billiard. 

3. Geometric character of the phase space structure, taking into account billiard border prop- 
erties (its diagonal, lacunas, discriminants). Universal character of symmetrical normal forms, 
describing the dynamics and local bifurcations in the neighbourhood of cycles and special zones of 
the billiard. 

4. Individualization of invariant distributions defined by the form of the billiard. The reduction 
of the measure (one-point factorization) without loss of statistic description. The division of quick 
and slow variables is not necessary. 

5. The unity of dynamics, kinetics and geometry of the billiard. The solution of direct and 
indirect problems on the restoration of involution, invariant measure and the form of the border. 
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Symmetric approach allows direct generalization on the multiply connected, multi-dimensional 
and other cases of different billiard border topology. The peculiarities named preserve here their 
key role. 
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Figure 

The geometry of a billiard and the schematic form of typical billiard symmetric phase space. 
Elliptic zones of regular motion R, a chaotic region C, the diagonal A, lacunas L and discriminants 
D. 
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